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We have developed a new boot-strap method for solving a class of interacting one-dimensional chiral 
fermions. The conventional model for interacting right-moving electrons with spin has an SO (4) 
symmetry, and can be written as four interacting Majorana fermions, each with the same velocity. 
We have found a method for solving some cases when the velocities of these Majorana fermions are 
no longer equal. We demonstrate in some detail the remarkable result that corrections to the non 
skeleton self-energy identically vanish for these models, and this enables us to solve them exactly. 
For the cases where the model can be solved by bosonization, our method can be explicitly checked. 
However, we are also able to solve some new cases where the excitation spectrum differs qualitatively 
from a Luttinger liquid. 

Of particular interest, is the so-called SO (3) model, where a triplet of Majorana fermions moving 
at one velocity, interact with a single Majorana fermion moving at another velocity. We show using 
our method, that a sharp bound (or anti-bound) state splits off from the original Luttinger liquid 
continuum, cutting off the X-ray singularity to form a broad incoherent excitation with a lifetime 
that grows linearly with frequency. 
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I. INTRODUCTION 

The anomalous normal state behavior discovered in 
the cuprate superconductors has stimulated enormous in- 
terest in the possibility of new kinds of electronic fluid 
that might provide an alternative to Fermi liquid be- 
havior. The classic model for non- Fermi liquid behavior 
is provided by the one-dimensional electron gas, where 
the generic fixed point behavior is a Luttinger Liquid. EJ 
Thanks to a wide array of non-perturbative techniques, 
there is a rather solid understanding of the Non-Fermi 
Liquid properties in such Id systems. Motivated by an 
early suggestion of Anderson,el many authors have at- 
tempted to generalize the Luttinger liquid concept to 
higher dimensions .ETu 

The Luttinger Liquid in Id is truly special in that it 
has no quasi-particle poles but a branch cut singular- 
ity; its correlation functions are scale-invariant, with an 
associated beta function that is zero to all orders in per- 
turbation theoryu for a wide range in the coupling: 

/3(g) = 0. 

That the beta function is zero is not in itself special to the 
Luttinger Liquid. For example, in the absence of nesting, 
or a Cooper instability, the beta function associated with 
Landau's Fermi Liquid fixed point is also zero for the 
forward scattering channelM 

The profound differences between the Luttinger liquid 
and Landau Fermi liquid fixed points originate in the spe- 
cial kinematics of one dimension. In Id, the Fermi surface 
consists of just two points ±fc/ where the electrons inter- 
act very strongly, and asymptotically near these Fermi 
points, energy and momentum conservation impose a sin- 
gle constraint on scattering processes, giving rise to a 



qualitative enhancement in scattering phase space. This 
causes the electron to lose its eigen-state status to the col- 
lective spin and charge density bosonic modes. Luttinger 
liquid behavior requires the absence of umklapp interac- 
tions, and in this case, left- and right-moving particles 
are separately conserved. The spin and charge current 
densities of the right (or left) moving particles are then 
simply proportional to the corresponding spin and charge 
densities: 
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so that the continuity equation assumes a special form 

(d r - iv s ,cd x )p R c = 



As noted long ago by Dzyaloshinskii and Larkintm, these 
conservation laws lead to the vanishing of the iV-point 
connected current correlation functions for N > 2 ( "Loop 
Cancellation Theorem", see Section IV), which leads to 
a Gaussian theory for the spin and charge bosons in the 
Tomonaga Luttinger model, and also for the low energy 
effective theory of the Hubbard model in Id. 

Unfortunately, the special kinematics of one dimension 
do not survive in higher dimensions, and largely for this 
reason, attempts to generalize the Luttinger Liquid to 
d > 2 with strictly local interactions have been unsuccess- 
ful. In one dimension, energy and momentum conserva- 
tion impose a single constraint on the forward scattering 
processes, whereas in higher dimensions, they impose in- 
dependent constraints on the scattering processes. These 
additional constraints eliminate many of the potentially 
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dangerous singularities present in one dimensional scat- 
tering processes, stabilizing the Eermi liquid in two or 
higher dimensions BU Lin et. alia arrived at the same 
conclusion, making the passage from one, to two dimen- 
sions by coupling N Hubbard chains together and taking 
the limit N — > oo. Whilst it is possible to circumvent the 
Fermi liquid in two dimensions by intrducing long-range 
or singular interaction aljE3 'BI, a route to non-Fermi liq- 
uid behavior in 2d that involves strictly local interactions 
has not yet been found. 

An alternative approach has however been advocated 
by Anderson,li3 who noted that higher dimensional non- 
Fermi liquid behavior might derive from the formation of 
bound, or anti-bound states above and below the single- 
particle continuum. Such bound-states play an impor- 
tant role in the formation of the one dimensional Lut- 
tinger liquid, where they give rise a finite scattering phase 
shift at the Fermi energy, driving the formation of X-ray 
singularities in the spinon-holon continuum. 

In this paper, we are motivated by this discussion to 
examine whether such singularities are robust against the 
removal of some of the special kinematic symmetries of 
one dimension. By modifying the Id kinematics, we show 
that it is possible to actually split-off bound-states from 
the spinon-holon continuum giving rise to a new type of 
one-dimensional non-Fermi liquid that does not rely on 
the special Id symmetries mentioned above. The key to 
our idea is as follows. The electron fluid on the Fermi sur- 
face is made up of spin- up and down electrons and holes. 
Borrowing from the Dirac equation, we can rewrite the 
electrons and holes as charge-conjugation eigen-states: 

C T = -L(*« _ ^(2)), C[ = --^(*< 3 > + i*<°>). 

where the \E , ' a - ) (a^ (0,1,2,3) ) represent four chiral 
Majorana fermionso such that ^f^(x) = ^ a ^(x). In- 
stead of changing the interaction, we modify the scatter- 
ing kinematics by making one of the Majorana fermions 
to have a different velocity to the others. In the classic 
Tomonaga Luttinger model, all four Majorana fermions 
have the same velocity, (exhibiting the full 5*0(4) sym- 
metry) and this leads to the special Id kinematics men- 
tioned. But in our model (with the reduced 5*0(3) sym- 
metry) , lifting the velocity degeneracy causes the energy 
and momentum conservation to be distinct constraints 
in scattering phase space. We shall show that in this 
case, the reduced (relative to the Luttinger model) scat- 
tering cuts off the X-ray catastrophe associated with the 
Luttinger liquid behavior. The "horn-like" feature in the 
spectral weight of the Luttinger liquid is then split into a 
sharp bound (or anti-bound) state that co-exists with an 
incoherent spin-charge decoupled continuum. We sum- 
marize these results in Fig.[j]. 

While the main motivation of our model has been to 
find new fixed point behaviour in Id, our model (0) also 
has physical relevance to two recent work: 

• The transport phenomenology of the cupratesEl 




FIG. 1. Schematic diagram showing the evolution of the 
spectral weight as we introduce velocity difference to the 
fermions. Inset indicates the bare spectral function, without 
interactions. 

suggests that electrons near the Fermi surface 
might divide up into two Majorana modes with 
different scattering rates and dispersion. To date, 
this kind of behavior has only been realized in 
impurity modeisL3 and their infinite dimensional 
generalization!!-] We shall show that by breaking 
the velocity degeneracy of the original chiral Lut- 
tinger model, we obtain a one-dimensional realiza- 
tion of this behavior: a sharp Majorana mode inti- 
mately co-existing with an incoherent continuum of 
excitations, reminiscent of the higher dimensional 
phenomenology. 

• Frahm et. al.@ have recently proposed that the 
low energy effective Hamiltonian of an integrable 
spin-1 Heisenberg chain doped with mobile spin- 
1/2 holes is given by (||), with one Majorana 
fermion describing a slow moving excitation 
coming from the dopant, interacting with three 
rapidly-moving Majoxana fermions that describe 
the spin-1 excitationsE3 of the spin-chain. (See Dis- 
cussion.) Such doped spin-chain models may be 
relevant to certain, experimental systems such as 
Y 2 _ x Ca x BaNi0 5 E] 

Whereas the, 50(4) model can be treated by 
bosonization,EI'E!J by changing the velocity of a single 
Majorana fermion we introduce a non-linear term into 
the bosonized Hamiltonian that preclude a separation in 
terms of Gaussian spin and charge bosons. (See Discus- 
sion.) 

To tackle this SO (3) model, we have developed a 
new fermionic bootstrap method, that has its basis the 
diagrammatic approach of Dzyaloshinskii and Larkin 
(1974)tl Their method depends crucially on the existence 
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of conserved currents to eliminate large sets of diagrams, 
leading to a closed set of equations that can be solved an- 
alytically for the Green function. On first glance, the re- 
duced number of conserved currents in the £0(3) model 
(compared to the SO (A) model) causes the Dzyaloshin- 
skii and Larkin method to be inapplicable, because one 
has to deal with non-conserved current vertices that in- 
volve the singlet Majorana fermion of different velocity. 
We have found however, that by dealing directly with 
fermionic propagators and the four-legs fermionic vertex, 
bypassing the intermediate currents, there are enough 
conservation laws after all to eliminate all vertex correc- 
tions to the skeleton self-energy (Fig. 2), allowing us to 
write down a compact set of coupled equations involving 
only the fully renormalized skeleton self-energy and the 
exact Green function of the theory. 

The plan of the paper is as follows. In Section II, we 
define the class of models of interest here. In Section III, 
we describe out, modification of the classic Dzyaloshin- 
skii and LarkinQ diagrammatic method for solving one- 
dimensional fermionic systems, to deal with our case 
where not all the velocities are the same. In Section 
IV, we take advantage of the purely chiral nature of our 
model (|^) to write down a scaling form to simplifying con- 
siderably the bootstrap equations derived in Section III. 
In Section V, we derive asymptotic solutions for frequen- 
cies near the spectral weight singularities, and demon- 
strate our results with numerical solutions. In Section 
VI, we discuss the nature of this new fixed point. Some 
of the results have appeared in a brief form in Ho and 
ColemanEl 



II. MODEL 

The class of model we study here is: 

f 3 
H= dx^-iJ2 v a^ {a) (x)d x ^ (a \x) 

+g^\x)^ 1 \x)^ 2 \x)^ 3 \x)Y (1) 

where "I^"' are real (Majorana) fermions such that 
^ a \x) = * (a)t (a;). The fermions are chiral (right- 
movers, say): this is a crucial property that ensures the 
system stays gapless, and allows for exact solutions in a 
number of cases. 

In the special case where all velocities are the same, 
this model has an 50(4) symmetry, where the four Majo- 
rana modes can be associated with the spin up and down, 
electron and hole excitations of the Fermi surface. To see 
this, write c T = ^(* (1) -M™), c x = -^(* (3) 
where c Q are the usual Dirac fermions, and the 50(4) 
model is just the conventional one-branch spin-1/2 Lut- 
tinger model: 




H. 



SO(4) 



j c er( x )2 v d x C a {x) + H.C. 



FIG. 2. Renormalized "skeleton self-energy" (SSE), where 
double lines represent full propagators. 



- g (4(aOc T (s) - 1/2) (c\(x)c L (x) - 1/2) }. (2) 

This 50 (4) modeLcan be shown by bosonization to be a 
Luttinger Liquid t3 

We shall mostly focus on the 50(3) model where v± = 
v 2 = v 3 = v ^ v : 

H = J dx^H kln (x) + .g*(°)(x)1'( 1 )(a;)* (2) (a;)* (3) (a;)}, 
3 

H km (x) = -tvY / ^ ia) (x)d x ^ a \x)-iv 'f (0) (x)d x ^ \x) 

a=l 



Note that this model reduces to the single-impurity 
model of Coleman et. al .113 when the mode ty(°> is 
made to localize at the impurity site, and Ho and Cole- 
man have studied the same lattice 50(3) model in high 
dimensions .EJ We will show that, by making the velocity 
of one Majorana fermion different, the scattering phase 
space decreases drastically, leading to this singlet split- 
ting off from the Luttinger continuum to form a sharp 
bound-state/anti-bound-state. Thus this is a system that 
has two qualitatively distinct relaxation rates, a dramatic 
departure from the Luttinger Liquid scenario. 

The 50(2) x 50(2) model where v = v x ^ v 2 = v 3 is 
also solvable by bosonization, and interestingly, our boot- 
strap method also works here. (See the sections Results 
and Discussion.) 

Finally, we shall also briefly look at the 50(2) model 
where vq ^ v\ ^ v-2 = v 3 . While we do not know if 
our method works here, we expect that due to the sep- 
arate energy and momentum conservation, there is still 
a restriction of scattering phase space, and the theme of 
split-off sharp bound/anti-bound-states continues. 

Note that the number of degrees of freedom and the 
interaction are the same in all the cases; the variety of 
behavior seen is due solely to changes in the scattering 
phase space, when the velocities of the fermions are made 
to be different. 



III. METHOD:PHILOSOPHY 

Our approach is based on the observation that for the 
50(4) and 50(3) models (and possibly others too), the 
renormalized skeleton self-energy, (SSE) containing full 
propagators, but no vertex corrections, (Fig. 2) is exact, 
so that 
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E (z, r) = g 2 G b {x, t)G c {x, r)G d (x, r), (4) 

where the G a are the exact, interacting Greens func- 
tions and {a, b, c, d} is a cyclic permutation of {0, 1, 2, 3}. 
These equations close with the usual relations: 

E (fc,w) = (^-w Q fc)-G a (fc,w)- 1 . (a = 0,1, 2, 3) (5) 

Equations (^,^) together define a boot-strap method to 
solve the problem. 

To show that there are no vertex corrections to the 
renormalized skeleton self-energy, we first review and 
then extend Dzyaloshinskii and Larkin's method. Pro- 
vided that we have a minimal 5*0(3) symmetry, then 
the three current densities j a (x) = — ie abc ^i^ {x)^!^ {x) 
(a,b £ (1,2,3)) are cor Wp|Wd classically. Following 
Dzyaloshinskii and LarkinntJeJ, since charge and current 
are the same in a chiral model, the continuity equation 
guarantees that the N-point connected current-current 
correlation functions vanish for N > 2: (xj = (2^,7$)) 

(j o (x 1 )i°(x 2 )...j«(x w )) c =0 ) (N>2) (6) 

For the non-interacting system, this result leads to the 
"loop cancellation theorem": for the amplitude associ- 
ated with a closed fermion loop with N > 2 conserved 
current insertions, the sum over all possible permutations 
of {xi} of the current operators must give zercotauJ. In 
Appendix A, for illustration, we give a derivation for the 
N = 4 case and also for odd N. Dzyaloshinskii and 
Larkin used this cancellation to eliminate all diagrams 
that contain such closed loops, considerably simplifying 
the vertex function and polarization bubbles. 

We use the loop cancellation theorem in a new way, 
to show that the vertex corrections to the skeleton 
self-energy (SSE) (Fig.2\-||ideptically vanish. Unlike 
Dzyaloshinskii and LarkirfltJ'Ea, we discard the interme- 
diate currents and the associated current vertices, and 
deal only with fermionic propagators and the four-leg in- 
teraction vertex. The Loop Cancellation Theorem is the 
same. This method has the advantage that it is more 
compact (only the self-energy and the Green functions 
are involved), and treats all propagators in a symmetric 
manner. To illustrate the idea, consider the self-energy of 
the singlet Majorana mode in the SO(3) model. Figs.|| 
list all such diagrams at order g . The Feynman dia- 
grams contributing to the skeleton self-energy are con- 
structed by combining loops with two insertions. This 
is clearly true for the second order diagram, and we il- 
lustrate this using the first non-trivial order: the fourth 
order diagram in Fig.[|(A), but it holds to all orders in 
perturbation theory. Non-skeleton contributions to the 
self energy involve diagrams with loops containing more 
than two current insertions. In these diagrams, the sum 
over all permutations of the current insertions into the 
loops is automatically zero, as illustrated to order g 4 in 
Fig.||(B). A convenient way to represent these diagrams 
is to split each diagram into a backbone which is the 



same in all three diagrams, and the 4-insertions loop. In- 
serting the 4 vertices of the 4-loop in various ways into 
the 4 vertices of the backbone gives the three diagrams 
in Fig.||(B). Note that this method of generating the di- 
agrams give rise to the correct degeneracy for each of the 
diagram types (i,ii,iii). 

To generalize these results to higher order graphs, it 
is more convenient to look at the set of diagrams for the 
free-energy. Cutting a ty( ' line gives back the singlet 
self-energy So- We first note that only even orders in 
g occur in the free-energy expansion, because the bare 
Majorana propagators are diagonal in the Majorana fla- 
vor index. Next, there is always a closed loop with n 
propagators (not necessary of the same type) in any of 
the free-energy diagrams of order g n . Otherwise, im- 
proper and/or disconnected self-energy diagrams would 
be generated. Then, at order g 6 for example, we have the 
following classes of diagrams listed in Fig^ that might 
generate non-SSE diagrams. 

The Loop Cancellation Theorem applies to each case 
where there is a closed loop with more than 2 propaga- 
tors of the same kind. Thus case (III) is the only one 
left. Yet, case (III) generates either SSE diagrams, im- 
proper self-energy diagrams (where cutting one of the 
lines lead to 2 disconnected parts), or else, diagrams that 
have already been counted in the other cases. The last 
observation follows from the fact one can always find a 
closed 6-loop or 4-loop buried in the diagram. Hence, all 
potential non-SSE generating diagrams disappear! One 
can clearly generalize the same reasoning to higher order 
diagrams. We only need to check that this method deals 
with the combinatoric factors correctly, i.e, all the degen- 
eracies of the diagrams are such that there are no non- 
SSE diagrams left over. Here, we appeal to the fact that 
in the 5*0(4) model, there must also be the correct loop- 
cancellations, because our method gives the same exact 
answer as Dzyaloshinskii and Larkin's method. Even 
though we have drawn the diagrams treating the triplet 
lines as identical, these triplet lines actually must carry 
a Majorana flavor index, and to generate all possible di- 
agrams whether distinct under SO (3) or not, we must 
draw all possible diagrams with proper indexing of each 
of the lines. Listing all diagrams this way is indepen- 
dent of which symmetry we are dealing with, and conse- 
quently, combinators factors will automatically be taken 
care of in doing loop cancellation with these Majorana 
indices on the propagator lines. In particular, the sym- 
metry or combinatoric factors for each diagram must be 
just right to allow loop cancellation to work in the 50(4) 
case, and hence for the SO (3) case too. 

Thus we can show that the vertex corrections to the 
self-energy So of the singlet Majorana fermion cancel 
to all orders, leaving the fully renormalized SSE as the 
only remaining contribution. Intriguingly, this argument 
fails for the SO (2) x 50(2) model, because each ver- 
tex has two "fast" legs and two "slow" legs, unlike in 
the 50(3) case where there are only one of the singlet 
leg. Thus, for example, the non-SSE diagrams in Fig.|| 
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(A) 




(i) + (ii) + (iii) = (Loop Cancellation Theorem) 

FIG. 3. (A) Illustrating how the only non-vanishing sin- 
glet self-energy at order g 4 is constructed by combining a 
propagator back-bone with loops containing two vertex inser- 
tions. Dotted lines indicate bare propagator for the singlet 
Majorana fermion ty^. Full lines indicate bare propagator 
for the triplet Majorana fermions. (B) Illustrating how the 
non-skeleton self-energy at order g 4 is constructed by combin- 
ing a propagator back-bone with loops containing four vertex 
insertions. 
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FIG. 4. List of all classes of free-energy diagrams that gen- 
erate non-SSE diagrams at order g 6 . 
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FIG. 5. Examples of non-skeleton self-energy diagrams in 
the 50(2) x SO(2) model where the Loop cancellation The- 
orem does not apply. 
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Ga(r/ix), 



(9) 



This form is motivated by the observation that chirality 
prevents space from acquiring an anomalous dimension. 
Under a Fourier Transform, this scaling form is self dual, 



^T—Ga{r/ix) —G a (k/iuj), 
Zmx lui 



(10) 



where the same function Ga appears on both sides. In- 
serting Eqnjl^ into ^| and Fourier transforming, 



1 



2nix du' 



[l-VaU-l/Ga{u)] u=T/tx . (11) 



Since the bare Green function scaling form is 1/Q®(u) = 
1 — v a u, it does not contribute to the self-energy. Com- 
bining Eqns. U and O, 



do not have a closed loop of only one kind of propaga- 
tor to allow loop cancellation to apply. Yet, these dia- 
grams cannot contribute to the exact self-energy either, 
because the SO{2) x 50(2) model can be solved exactly 
by bosonization, or by a slight extension of Dzyaloshin- 
skii and Larkin's method, and these results agree exactly 
with our bootstrap method (See Discussion). There must 
then be more cancellation than due to just the Loop Can- 
cellation Theorem in its current form. 

To complete our proof, we need to show that the triplet 
Majorana self-energy is also given by the skeleton dia- 
gram. We use the full Kadanoff-Baym Free energy Func- 
tional: 



F[G] = -r{Trln[G _1 ] + Tr[£G]} + Y[G], 



(7) 



where Y[G] is the sum of all skeleton diagrams.E-3 Now, 
by construction, SF[G]/5G a = generates the equations 
for the self-energies, and in particular, 6F[G]/5Gq must 
generate the skeleton self-energy Eo. This requires that 
the the Kadanoff-Baym Free energy Functional truncates 
at the leading skeleton diagram: 



F = -TjTrlntG- 1 ] + Tr[£G]} 




(8) 



Finally, by differentiating the Free-energy functional with 
respect to the exact Greens functions 01,2,3 of the triplet 
Majorana fermions, each triplet self-energy is also given 
by the corresponding skeleton self-energy. 



IV. METHOD-DETAILS 



^2 ISaiu)}- 1 = -(g/27r) 2 G b (u)Gc(u)Gd(u) (12) 

where {a,b,c,d} are cyclic permutations of {0,1,2,3}. 
The boundary conditions are: 



0.(0) = 1, G' a (0) = 



(13) 



derived from the physical requirement that at high fre- 
quencies, the fermions are free particles, moving with the 
bare velocity v a . Equations [l^ and |l3| are the scaling form 
version of our bootstrap method equations fi|, ^[ Notice 
that the differential equation [12], like Eqn.[I|7is indepen- 
dent of the sign of the coupling g. Also, ( |l2| ) has no 
information on which model of the class (|j) it refers to, 
the symmetry of the model (ie. the velocities) only come 
in through the boundary conditions (|l3|). 

For the 50(4) model, where G a {u) = G{u) (a = 
0, . . . , 3), Eqns. [l2| reduce to a single differential equa- 
tion: 



^[awr^-G^Tr) 2 ^)] 3 , 



(14) 



for which the solution satisfying the boundary conditions 
g(0) = 1, G'(0) = v is: ' 



G(x,t) 



1 



2nix 



[1 — v+r/ix] [1 — v-r/ix] (15) 



where v± — v ± (g/2n) and v is the bare velocity. Iden- 
tical results are obtained by bosonizatiorO, where i>+ 
and V- are in fact the velocity of the spin-boson and 
the charge-boson. Thus this confirms that the skeleton 
self-energy is exact for the Luttinger model. 



We now apply this result, using the limiting case of 
the 50(4) model to check the validity of our results. Our 
equations are dramatically simplified by seeking solutions 
to (0) which satisfy a scaling form 



V. RESULTS 

In the 50(4) model, the electron spectral weight dis- 
plays two classic X-ray singularities associated with the 
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FIG. 6. Spectral weight of SO (4) model 




FIG. 7. Spectral weight of the SO(3) model. For clarity, we 
have shifted up the curves for various momenta by 8 units. 



decay of-the electron into a spinon and holon continuum. 
(Fig.^).EJ We now show that if Aw — v — vo is finite, 
one of these X-ray edge singularities is completely elimi- 
nated. If wo < w, we find that low velocity "horn", origi- 
nally with velocity w_ , develops a sharp bound-state pole 
in the singlet channel, and a broad incoherent excitation 
in the triplet channel with a lifetime growing linearly in 
energy. If vq > w, the high velocity "horn" splits off 
a singlet anti-bound-state and the triplet channel devel- 
ops a high- velocity incoherent excitation. (Fig.Q). The 
sharp bound-state in the singlet channel develops once a 
velocity difference is introduced, because energy and mo- 
mentum conservation now provide distinct constraints to 
scattering (unlike in the 50(4) model), leading to much 
less phase space for \I/( ) to decay into. 

To see this, we must analyze Eqn.[l2] for the 50(3) 
case: 
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FIG. 8. Quasi-particle weight Z of * (0) in the SO (3) 
model. 



d 2 



du 2 



g^ 1 = -(g/2n) 2 (G 3 ) 3 



(16) 



A very convenient way to discuss these equations is to 
map them onto a central force problem. If we write r = 
{Qs^Go 1 ), F = -{gg 3 /2ir) 2 {g ,g 3 ), then, f = F, where 
r = ie., u is like "time". By inspection, r x F = 
0, so the force is radial, thus the "angular momentum" , 
r x r = Aw is a constant. If we use polar co-ordinates, 
(C/" 1 ,^" 1 ) = r(cos 9, sin 9) the equations for the Green- 
function resemble the motion of a fictitious particle under 
the influence of an anisotropic central force: 



Av 2 



= -(g/2n) 2 - 
= Av. 



1 



' sm f 



(17) 



The velocity difference Aw — v — vq provides a repulsive 
centrifugal force. The boundary conditions (|l^) mean 
that the "particle" starts out at r(0) = y/2, 9(0) = tt/4, 
and with a slope change 9(0) — Aw/2. 

Without loss of generality, let Aw < 0. For Aw > 
simply replace w+ — > w_, g — * —g. When Aw = 0, the 
"particle" falls directly into the origin, and both g 3 and 
go diverge with X-ray singularities when the particle first 
hit the origin at "time" u = l/w+. Then the particle goes 
pure imaginary in both coordinates, which gives rise to 
the Luttinger continuum in the spectral weight, until the 
time u = 1/v- when the particle goes back to the origin, 
leading to the other x-ray singularities for both g 3 and 
go. From then on, the particle stays in the real plane. 
(Fig-i 

However, once Aw < is finite, 9(0) = Aw causes the 
orbit to miss the origin at u ~ V u ±, - Instead, 9^0 
at some finite "time" u = 1/wg (Fig.O), at which r = C 
and 9 = Av/C' 2 . For u ~ it follows that (r,9) = 

(C, 9 (u — 1/wg)) , from which we can read off the following 
asymptotics: 
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Gziu)- 1 -^ (18) 
g {u)- 1 ~{l-uvS)/Z, Z = Cv* /\Av\. (19) 

Thus the associated X-ray singularity in the spectral 
function for both G3 and Go is eliminated, replaced by an 
anti-bound state for the singlet Go with spectral weight 
Z, moving with velocity Vq , splitting off above the contin- 
uum. After this time, r is complex in both coordinates, 
until eventually at u — the particle passes through 

the origin, giving rise to the remaining x-ray singularity 
at it = 1/^3 in both G3 and Go. 

The quantity £ = (v — vo)/g plays the role of a cou- 
pling constant, and approximate analytic solutions are 
possible in the limiting cases of small and large £. For 
|At>| >> Jj^ interactions can be ignored, so — > vq, and 

Z — > 1~. For I Aw I << the "motion" of the fictitious 
particle emulates that of the 50(4) model until the an- 
gle 9 approaches zero. We may estimate Vq and C by 
integrating Eqnjll] with the approximation r(u) w r(u), 
where f = [2(1 — v+u)(l — v-u)} 1 ^ 2 is the 50(4) solution: 



du 



-du 



7T 

T 



C « r (IK). 



(20) 



After doing the integral, this estimate gives (for |Av| << 
M/27T): 



v Q = v + H — exp — 

7T 



.9 



Z = 



V2 5 



7tA« 



exp - 



2Av 
3 



4Aw 



(21) 
(22) 



indicating that the formation of the sharp anti-bound- 
state is non-perturbative in the velocity difference. 
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FIG. 10. Spectral weight of the 50(2) x SO(2) model. 



To illustrate these results further, we have carried out 
numerical solutions of the differential equations ([l6|) for 
intermediate values of the coupling constant £, using a 
standard adaptive integration routineHj Results are sum- 
marized in Figsj^,^. 

While we have not established the validity of our 
method to models of lower symmetry (but see Discus- 
sion), we believe that the method captures the essence 
of the kinematic constraints imposed by energy and mo- 
mentum conservation, at least for weak coupling. Thus, 
we have also performed numerical calculations for the 
SO(2) x SO(2) and the 50(2) models. 

For the SO (2) x 50(2) model, the pair $(°), tfC 1 ) with 
the same bare velocity can combine together to form a 
boson, and similarly for $( 2 ) ) ^>( 3 ). This leads back to a 
Luttingcr liquid form, but with asymmetric power law 
singularities at the renormalized velocities t>+ and V- 
(Fig.|To|) . (Also see Eqn.29 later in Discussion, for the 
exact analytical solution for this model.) 

As we progress to the SO(2) case, when vq < v\ — 
v 2 < W3, we see a sharp pole for the fermion which has 
the extremal velocity different to all the others, while the 
Luttinger continuum turns into wide peaks linear in en- 
ergy for the fermion(s) with intermediate velocities, see 
Fig.|ll|. This illustrates once more our contention that 
making one Majorana degree of freedom to have a dif- 
ferent (extremal) velocity causes drastic collapse of the 
scattering phase space for this fermion. 



VI. DISCUSSION AND CONCLUSION 

A. The Id Majorana SO (3) Model 

In summary, we have demonstrated that by breaking 
the velocity degeneracy of a system of interacting chi- 
ral fermions we restrict the scattering phase space in a 
way which causes a sharp bound or anti-bound state to 
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FIG. 11. Spectral weight of the SO (2) model. For clarity, 
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split off from the spin-charge continuum, leading to a 
system with two qualitatively distinct spectral peaks and 
scattering rates. This is a significant departure from the 
Luttinger liquid scenario and demonstrates a new class 
of one-dimensional fixed point behavior. 

This new fixed point exhibits properties in common 
with both Luttinger and Fermi liquids, and is per- 
haps closest in character to the Marginal Fermi liq- 
uid phenomenology introduced in the context of cuprate 
metals.E3 Like the Fermi liquid, there is a sharp quasi- 
particle bound-state, but this co-exists with a Luttinger 
liquid-like continuum which is bounded by two extremal 
velocities. 

As mentioned, the 50(3) model cannot be solved by 
conventional bosonization, forcing us to introduce this 
new bootstrap method. Two immediate questions arise: 

• the nature of the 50(3) fixed point, and 

• the range of validity of the bootstrap method. 

In the SO (A) model, the fermionic spectral weight has 
X-ray singularities at the velocities v + , V- (see Methods- 
Details). By bosonization, the model can be mapped 
onto a theory of free bosons (the spin-boson and charge- 
boson) moving at v + ,t>_, where for g > 0, v sp i n — V-, 
Vcharge = v +, and for g < 0, the role of are 
swapped round. This is a direct consequencej-|Of sepa- 
rate charge and spin conservation in the model.u We can 
demonstrate this in the Majorana fermionic representa- 
tion. The classically conserved densities are: 



J 01 (x) = -i^ i0 \x)^(x), 
J 23 {x) = -WW(x)¥ 3 \x). 



(23) 



(By the 50(4) symmetry, we can define other combina- 
tions also.) Using the commutation relations listed in 



Appendix B, we get the equations of motion: 



{-d T - vq)J Q1 (q) = —qJ2s{q), 
{-d T - vq)J 23 {q) = y-qJoiiq)- 



(24) 



The right hand side of the equations is not zero (as 
would be expected for conserved currents) because of the 
anomalous commutator (Appendix B): 

[Joi(p), Joi(g)] = [J 23 (p), J2 3 (q)} = P 6(p + q), (25) 

which is the SU (2) level 2 Kac-Moody algebra anomaly.0 
Fortunately, by diagonalizing the system (p4]). the lin- 
ear combinations J-(q) = Jqi(q) — J2 3 (q) and J+(q) = 
Joi(q) + ^23(9) do satisfy the continuity equation: 



(-d T - v+q)J+(q) = 0, 
{-d T - v-q)j_(q) = 



(26) 



where v± is as before in Eqn.(|l5|), indicating that these 
new densities rJ± are proportional to the spin-boson and 
charge-boson.Ell This then leads to sharp poles in the 
charge and spin susceptibilities. 

For the 50(3) model, using the same definitions (|22 
we find: 



{-d T - v q)J i(q) = ^-qJ23(q) + («3 - v )K 01 (q), 
{-d T - v 3 q)J 23 {q) = — qJoi(q). 

Z7T 



(27) 



where ^01(9) = -iJ2 k k ^ (0) (q ~ fc)* (1) (fc)- This extra 
term came from the commutator of J01 and the kinetic 
energy, and causes the set ( p7p not to close, and bosoniza- 
tion in terms of free spin and charge-bosons (or any lin- 
ear combinations) is impossible. In short, because of the 
anomaly, the classically conserved 5*0(3) density J 23 is 
admixed with the classically non-conserved J01, leading 
to the loss of a sharp pole for the susceptibility corre- 
sponding to J 23 . Thus it is unlikely that the model can 
be mapped to a model of free bosons, which makes it very 
different to the conformally invariant fixed points of the 
SO(4) model and the 50(2) x 50(2) model (see below). 
Also the presence of a sharp pole in the fermionic spec- 
tral weight indicates that there is at least one fermionic 
degree of freedom in -the diagonalized system. 

Now, Frahm at.al.E3 has conjectured that this 50(3) 
model is the low-energy effective theory of an integrable 
model of a spin-1 chain doped with spin-1/2 mobile holes. 
Using Thermodynamic Bethe Ansatz, they have shown 
that the spin and charge sectors of the doped holes be- 
come decoupled at low temperatures, and have calculated 
the low temperature free energy of the spin contribution 
to be: 



7rT 2 

6v 
ttT 2 

6v 



1 3A 

2 ~ 4^ 

3 3A 
2 + 4^ 



In A 



In A + 



(28) 
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where A > is a constant that depends on the doping 
only. Withr^ = (undoped case), the first term has been 
interpretedliS as coming from a single Majorana fermion 
of velocity vq, while the second term come from a triplet 
of massless Majorana fermions with velocity v that rep- 
resent the SU(2\-leve\ 2 WZNW model, which has been 
shown by Afflecktj to be the low energy effective theory of 
the gapless integrable spin-1 chain. Naively, one expects 
that a system of fermions with two velocities cannot be 
conformally invariant, unless the two species do not inter- 
act with each other and thus form two decoupled sectors 
that are individually conformally invariant S3 If the free 
energy ( p8| ) is indeed the 50(3) model low temperature 
free energy, the form of the free energy suggests that the 
SO (3) model is again conformally invariant asymptoti- 
cally, and hence the diagonalized (presumably fermionic) 
basis consists of two decoupled sectors. Note that the 
diagonalized basis is unlikely to coincide with the bare 
Majorana fermions of Eqn.(|3|)-— because the coupling is 
marginal, at least up to 0(g^)t3 

As for the range of validity of the bootstrap method we 
introduced to solve the model, we note that if we change 
two Majorana velocities at the same time, so that vq = v\ 
and x>2 = v 3 , we would have reduced tho-symmetry still 
further, to an 50(2) x 50(2) symmetry.^!! We can solve 



the differential equations (12) with the results 



G 3 (x,t) = 



2nix 
1 



V + T 

ix 



1 - 



V-T 



V± = - [Vq + V 3 ± [(v 3 - Vq) + (g/Tr)' 



(29) 



and 7 = ±(v 3 - v ){(v 3 - v ) 2 + {g/ir) 2 } a . Interest- 
ingly, this model can be bosonized to a model of free 
bosons, and the bosonization result agrees exactly with 
(^9|). This is surprising because as far as we can see, the 
closed-loop cancellation is not sufficient in the case of the 
50(2) x 50(2) model to cancel all vertex corrections. 
This suggests that a more general cancellation principle 
is at work, and that the range of validity of our solution 
may even extend to models with a still smaller, 50(2) 
symmetry. To date, we have not been able to prove this 
result. 

We also wish to point out that our differential ver- 
sion ( |l2| ) of the bootstrap equations (QJ^) are of such a 
simple form only because we have a purely chiral sys- 
tem. If we allow Left- and Right- ppvers to interact, 
the scaling form (||) no longer appliesEil, and we have not 
found a different scaling form that allows similar simpli- 
fications. However, we expect the bootstrap method to 
work still, as long as there are separate conservation of 
Left- and Right- currents. This is true at least for the 
50(4) model, because Dzyaloshinskii and LarkinO have 
shown that their method works also for such systems, 
and our method is a generalization of theirs. 



= v 2 = ^ 

SO(4) 



v i = ^> = ¥ 
vo f= M 

SO(3) 



1 d 



Luttinger Liquid 



spinon-holon 
continuum 



"Majorana Liquid 

J 



(anti-) 

bound-states 



2d 



Fermi Liquid 



sharp 

quasi-particles 



FIG. 12. 



B. Broader Issues: higher dimensions? 

Our work raises the question whether this kind of 
non-Fermi Liquid behavior might survive in dimensions 
higher than one. In higher dimensions energy conser- 
vation and momentum conservation are distinct con- 
straints on scattering phase space, and the Luttinger 
liquid rcvcrtfL.to a Fermi liquid, at least for short-range 
interactionsoB. By contrast, the SO (3) model cannot be 
solved by bosonization, and its unusual properties have 
reduced reliance on the special kinematics in Id. Thus, 
this kind of behavior might be more robust-, in higher 
dimensions. In fact, near infinite dimensionai-1, two life- 
times behavior persists in the 50(3) model, but here, 
the thermodynamics near zero temperature is that of a 
Fermi Liquid. The case of small, but finite dimensions is 
however, still open. 
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VII. APPENDIX A 

In this appendix, we prove by diagrammatic method, 
the Loop Cancellation Theorem for a loop with four cur- 
rent insertions. It is easiest to prove this in x,t space. 
(For a proof in momentum-frequency space, see Kopi- 
etz et.alrJ.) Let the four insertions be at x$ = (xi,Ti), 
i = 1, . . . , 4. Each leg of the loop is a free propagator: 



Gij — G(x^ , Ti , Xj , Tj ) 



1 



+ i{xi 



(30) 
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Denote by [1234] the loop where going clockwise starting 
from xi, we encounter successively xi,X2,X3,X4, ie., 



[1234] = G 43 G 32 G 2 iG 



14- 



(31) 



Without loss of generality, we can fix Xi and sum over 
permutations of the other three vertices. The Loop Can- 
cellation Theorem then says: 

[1234] + [1243] + [1342] + [1324] + [1423] + [1432] = (32) 

But for even number of propagators in a loop, going 
clockwise is the same as going anti-clockwise, hence, e.g. 
[1243] = [1342]. So, we only need to prove: 

[1234] + [1243] + [1324] = 0. (33) 

To do this, we need the important identity: 



GijGjk — Gik(Gij + Gjk), 



(34) 



which can be proven simply by substituting in Eqnj30[ 
Use this to rewrite the loops: 

[1234] = G14G43G32G21 = Gi3(Gi4 + G43)G32G21 
[1243] = G13G34G42G21 = Gi3(G34 + G42)G32G2l 

[1324] = G14G42G23G31 = G3i(Gi4 + G42)Gi2G23 (35) 

and it is clear that they do all cancel, since Gij = —Gji. 

^From this example, we can see that it is important 
for the cancellation of loops with even number of current 
insertions, that all the propagators be of the same type, 
to use the identity (134]). In our context, this means all 
the propagators are for fermions of the same velocity. 

For an odd number of insertions, the identity (34) is 
not needed, because time-reversal invariance guarantees 
the cancellation: a loop [lijkl . . . xyz] will be cancelled 
by the counter-clockwise partner [lzyx . . . Ikji], thanks to 
Gij = —Gji and a total of odd number of propagators. 
(This is the analogue pf Furry's theorem in QED, see eg. 
Peskin and Schroederc3.) 

We note in passing that this identity (Q) can also be 
used to prove the Ward Identity by diagrammatic meth- 
ods, order by orderB3 (For further information on the 
Ward Identity and how it is used for diagrammatic meth- 
ods for finding the exact Green Function in some one- 
dimensional systems, see Metzner et.alfl.) 



VIII. APPENDIX B 

Here we list some commutation relations used to de- 
rive the equations of motion of the various (classically) 
conserved densities. Start from the canonical anti- 
commutation relation for the Majorana fermions: 



(36) 



With the definitions (|23j), and with the SO (3) Hamilto- 
nian H = H Q + H int : 



f 3 
H = dx-i^2v¥ a) (x)d x ^ {a) (x) 

0=1 

-iv^ {0) {x)d x ¥°\x), 

Hi„t = -9 dxJ 01 (x)J 23 (x) (37) 



we can get straightforwardly: 

[J 01 (p),H ] = vpJ 01 (p) + (v- v )K Q x(p), 
[J 23 (p),H }=vpJ 23 (p). (38) 

We can recover 50(4) results by setting v = vq. Ordinar- 
ily, we would expect [Joi(p), Hi nt ] = 0, but this is spoilt 
by the SU(2) level 2 anomalous commutators: 

[Jbi(p),Jbi(g)] = [J r 23(p),J23(g)] = P S(p + q). (39) 

One can derive this by eg. a diagrammatic method, see 
Ch. 13 of the book by TsvelikB. This then leads to the 
only non-trivial commutation relations: 



[Joi(p),H lnt ] = ^pJ 23 (p), 

[J23(p),H lnt ] = —pJoiip)- 



(40) 
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